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Q Abstract 

Investigate the stochastic dynamic non- linear system with the Wiener 
and the Poisson perturbations: 



dx(t) = [P(t; x(t)) + Q(t; x(t)) • s(t; x(t))j dt+ 
+e (t; x(t) , dw(t) ,v(dt; d'y)) , 

< ^ where x G M n , n > 2; e(-) - is a strong random perturbation; w(£) - 

is the m-dimensional Wiener process; v[t; A7) - is the homogeneous 
on t non-centered Poisson measure. For such systems we construct 

<^ the program control s(t;x(£)) with probability one, which allows the 

CN system (1) to move on the given trajectory: itj(i;x(i)) = 0, i = 1, k. 

T 7"! Control program s(t; x(t)) is solution of the algebraic system of linear 

equations. Considered algorithm is based on the first integral theory 

^ for stochastic differential equations system. 1 

S3 

1 Preliminary definitions and results 

Suppose that a dynamic system is affected by random disturbance, which 
takes the form of Wiener and Poisson perturbations. In such situation we 
require that important system properties are kept invariant. We investigate 
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the program control construction problem with probability one for the similar 
systems. We use the first integral concept for the Ito's stochastic differential 
equation system (SDES) [1, 2, 3] and the automorphic function construction 
[4,5]. 

Let us consider the random process x(t) = x(t, u;), x : [0; +oo) — > 
which is the solution of the Cauchy problem for system of the Ito's generation 
stochastic differential equations (GSPE) 

m _ 

dxiit) = a,i(t; x(t))dt + bik(t] x(t))dw k (t) + / gi(t; x(£); ^)v{dt] d*y) 

A— 1 

R(<y) 



x(£) = x(£, x D ) = x , i = l,n, t>0, 
lt=o 

(2) 

where w(t)) is the m-dimensional Wiener process, z/(£; A7) which is homo- 
geneous on t non centered Poisson measure. This equation we present in the 
vector form 



dx(t) = A(t;x(t))dt + B(t;ic(t))dw(t) + J v(dt\d^) ■ G(t;x(t);j) 

B( 7 ) 

We suppose that the coefficients a(t; x), b(t; x), and g(t; x; 7) are bounded 

, . , .11. ! • dai(t:x) dbik(t;x) 
and continuous together with their own derivations — , — '— , 

OXj OXj 

— ' — i,j — 1, n respect to (t;x;7), and they are satisfied the condi- 

dXj 

tions of the theorem of existence and uniqueness of the Eq. (2). 

In the article [l] the concept of the first integral for Ito's GSDE system 
(without Poisson part) has been introduced. In the article [ ] the concept of 
the stochastic first integral for Ito's GSDE system with the centered measure 
has been given. Let's introduce the similar concept for the non-centered 
Poisson measure. 

Definition 1 Let random function u(t;x;oj) and solution x(t) of the Ito's 
GSDE system (2) be defined on the same probability space. The function 
u(t] x; (J) is called the stochastic first integral for Ito's GSDE system (2) with 
non centered Poisson measure, if condition 



u(t;x(t;x(0));w) = u(0;x(0);u 
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is held with probability equaled to 1 for the each solution x(£;x(0);u;) of the 
system (2). 

The random function u (t; x(£)) G C t f defined over the same probability 
space as the random process x(£) which is the solution for the system (2) 
is the first integral of the system (2) iff the function w(t;x(t)) satisfies the 
terms C) [6]: 

du(t \ x) 

1. &ij.(£;x) — — = 0, for all k = l.m (compensation of the Wiener's 

OXi 



2. 



perturbations); 

du(t;x) <9w(t;x)r 1 96j fc (t;x 



dt dxi 
of the time); 



Oj(i;x) — - bjk(t;x) — - — - — = (independence 



3. u(t] x) — u(t\ x + g(t; x; 7) j = for any 7 G -R(7) in the whole field of 
the process definition (compensation of the Poisson's jumps). 

Remark 1 // w;e analyze the concrete realization then the non random func- 
tion u(t;x) is the determinate first integral of the stochastic system. 

The concept of the stochastic first integral for the centered Poisson measure 
was introduced in [3]. The obtained conditions for its realization take the 
necessity of determining the intensive Poisson distribution density. In this 
research that condition is missing. Thus, it makes no difference what is the 
probability distribution of intensities of Poisson jumps. This case is very im- 
portant for constructing program controls [5]. In case a Wiener disturbances 
only, we construct the program control with probability one in article [7] on 
basis the research [2]. 

Theorem 1 [5] Let the determinate function u(t,x) be continuous together 
with their own derivatives with respect to all variables u(t, x) and the ran- 
dom function u(t, x; u) defined over the same probability space as the random 
process x(t) which is the solution for system 



dx(t) = A(t;x(t))dt + B(t;x(t))cZw(t) + J u(dt; dj) ■ G(t;x(t);i) 

R(j) ( 3 ) 

x(t) = x(£,x G ) = x , t> 0, 

t=o 
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where x G W 1 , n > 2; w(t) is an m- dimensional Wiener process, u(t; A7) - 
is homogeneous on t non centered Poisson measure. Assume the vectors e a; 
ex, e n consist an orthogonal basis in IR + x W 1 . If function u(t, x; u) is the 
first integral of the system (3), then the coefficients of an Eq. (3) and the 
function u(t, x) are related by conditions: 

n 

1. coefficients Bk(t;x) = x)e« (k = l,m), are columns of the 

i=l 

matrix B(t;x), that belongs to a set of functions 

( 



B k (t;x) G < 



q 0O (t; x) • det 



ei 




du (t; x) 


du(t; x) 


dxi 


Ox n 


/31 


fzn 


fnl 


fnn 



\ 1 



(4) 



where q 00 {t] x) is an arbitrary non-vanishing function; 
2. coefficient A(t; x) belongs to a set of the functions, defined by 



\ 1 n 

A(t;x) e\R(t;x) + -J2 



k=l 



<9x 



■BJ-) 



(5) 



where a column matrix x) components rj(t;x), i = l,n, are 
defined by: 

n 

C~ l (t\ x) ■ det x) = e + r^i; x)^; 

i=l 

C(t;x) an algebraic adjunct of the element e*o of a matrix if(t;x) 
and det C(t; x) 7^ 0; a matrix H(t; x) is defined as 



H(t;x) 



e 


ei 




du(t; x) 


<9w(t; x) 


du(t;x) 






dx n 




^31 






^n+1,1 


h 

• • " , n+l,n 



(6) 



and 



ggfc(^x) 
<9x 



zs a matrix of Jacobi for function B^it; x); 
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3. coefficient G(t;x;j) = gj(t; x; 7)^ belonging to Poisson 



measure, 



i=l 



is defined by the next representation G(t;x;7) = y(t;x;7) — x, where 
y(t; x; 7) is the solution of the differential equations system 



det 



c?7 

ei e 2 
aw(t;y(-;7)) dn(t;y(-;7)) 
<9yi dy 2 

<P3i{t;y (■■,!)) v ? 32(t;y(-;7)) 
^m(*;y(-;7)) ^n2(*;y(-;7)) 



du(t;y(-;i)) 
dy n 

Psnfay {■■,!)) 



(7) 



This solution satisfies the initial conditions: y(t; x; 7) 



x. 



7=0 



The arbitrary functions fj = /;.,•(£, x) ; hij = /^(t,x), ^ij = fij{t; y(s 7)) 
are defined by the equalities fij(t,x) = ^ - — hij(t,x.) = — ^— ^- — 

VyC^yC-^)) = ' and sefo of the functions j/jj, j/tjj, jy?, 

and the function u(t; x) together consist of the class of independent functions. 

Proof. We consider three steps for the proof. 

1. Let us use the first statement from the conditions £): 



1=1 



| du(t; x) 
dxi 



0. 



for all k = l,m. 



(8) 



<9it(£; x) 



e*i hold, then Eq. (8) 



If B fc (t; x) = 6«(t; x)ej and V x u(t; x) = ^ 

i=l i=l 

is an orthogonal property of vectors Bk(t;x.) and V x u(t;x). 

Taking into account a vector product definition in R n and their properties 
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we set a representation for columns Bk(t; x) of matrix B(-) — (^Bi(-), . . . ,B n 



B k {t;x) E < 



q 00 (t; x) • det 



/ 


ei 


■ ■ &n 




> 




du (t; x) 


du(t; x) 








dx\ 


dx n 








/31 


■ ■ f'in 




> 


V 


fnl 


■ ■ fnn 


J 





(4) 



where f\ = /j(t,x), i — 3, n, and functions x) = ^ - — - and u(t,x) 

together generate the collection of a independence ones. 

2. Later we use the second statement of the conditions £): 



dt 4-^ dxi 



8=1 



a»(*; x ) - ^L M^ x ) 



fc=i j=i 



0. (9) 



Let us suppose 

Q(t; x) = 1 + £ a j( t; x) - i £ £ £ x ) * ' ^ ' ? " " J 



t=i 



i=i fc=i j=i 



(9s i 



According to the scheme in article [2] we consider two vectors: 



□u(t; x) = ^^ X) e 



+ ^^x) e 



i=l 



and 



n m n 



<J(t;x) = e + ^a 4 (t;x)e; - | ^ 6j fc (t; x) dhlk ^' x ) e , 



i=l 



0Xj 



8=1 fc=l j = l 

Then the Eq. (9) means that the vectors Uu(t;x) and £^(t;x) are or- 
thogonal. Later we use the vector product and it's properties again, and we 
obtain: 



x e < 



det 



e 


e*i 


■ ■ &n 


<9w(t; x) 


<9w(t; x) 


du(t; x) 




dxi 


8x n 


^30 


hi 


hn 


hn+1,0 


hn+1,1 


• • hn+l,n 



>={detH}, (10) 
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where hi = /j(t,x), % = 3,n + l, and functions hij(t,x.) = — ^— - — and 

OXj 

u(t,xi) together generate the collection of a independence ones. Let us con- 
sider the vector 

n j n m n db (t~ X.) 

A(t;x) = e Q + ^ai(t;x)ei = Q(t;x) + - b jk (t;x) lk — ej. 



j=i «=i fc=i j=i 

As a coefficient of e G is equal to 1, then we get: 



n m n 



A (t ; x) G I C- 1 (t; x) • det H(t ; x) + - ^ ^ ^ 6, fc (i ; x) 



i=i fc=i j=i 



db ik (t;x.) 



where C(t;x) is an algebraical adjunct of element e of matrix H(t;x), 
detC(t;x) ^ 0. As far as the vector C _1 (i;x) • detif(t;x) we present in 
a form: 

n 

C _1 (£;x) • det if (i;x) = e + yVj(t; x)e^, 

i=i 

let us introduce the following vector x) that has components rj(i;x), 
i — l,n. 

We might the next representation: 



n m n 



i=i fc=i j=i 



dbik&x) 

dx,; 



= E 

fe=i 



fl£ fc (t;x) 

<9x 



S fc (*;x), (11) 



where 



<9x 



is a matrix of Jacobi for function x). It follows that 



a coefficient x) is a sum of matrices: 



9B fe (*;x) 



<9x 



£ fc (t;x) > , 



A(t;x)eii2(t;x) + ^^ 

L fe=i fc=i - 

3. According to the third statement in £) for all 7 G #(7) we have: 



w(t; x; w) — -u^t; x + x; 7); = 0. 



(12) 



It means that the function u(t; x; 00) is automorphic function under transla- 
tion by x with function G(t; x; 7). Let us set a condition for it. 
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According to [3] we have y(t; x; 7) = x + G(t; x; 7). For notational sim- 
plicity we take off a parameter u. Then we obtain 



w(t;x) = uft;y(t;x;7)j 
Ou(t] x) 

for all 7 G RH)- Hence, — — - — = and it holds: 

£77 

9u(*;y(*;x;7)J ™ <9u(t;y(t;x;7)J ^(t ;x;7 ) 



(13) 



$7 



i=l 



$7 



0. 



(14) 



n 5uft;y(-;7) 

Eq. p denotes that the vectors V y K.(i;y(-;7) ) = / > ~ ^ an( ^ 



i=i 



dyi 



^- ' ^ = ^- ' ^ e*j are orthogonal and they connected by the relation 

c>7 ^-^ oj 

i=i 



f>7 



G < 



det 



ei 

0w(t;y(-;7) 
dyi 



c n 

du[t;y(-;j) 
dy n 



• (15) 



There are (pi(t; y), i = 3, n such functions that (fij = ^ and u(t; y(-; 7) 

together generate the collection of a independence functions. 

By virtue of the fact that y(t; x; 7) = x + G(t; x; 7) then the Eq. (15) 
is the differential system, which has a function y(-;7) as unknown. Let 

<9y(-; 7) 

us expand the determinate (15) into the first row. Then we have — — — = 

07 

n 

Si(y (■; 7)) efj, where a is an arbitrary function, which independent from y. 

=1 

Such away we obtain the next differential system: 
f dy l {-- 1 ) 



a 



i=i 



c?7 



"5i(y(-;7)), 
a^n(y(-!7))- 



(16) 
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Suppose that y(i; x; 7; 9) is a solution of the Eq. (16), where 9 is a constant 
vector. By virtue of (13) holds for all t, x and for all 7, then is 

u(t;x) = u(*;y(t;x;7i;0)J = w^;x+G(t;x;7i; 0)J = u(t; x+G(t; x; 7 2 ; 1 



(17) 

As a special case for same 7 = 7 the Eq. (17) is determined by expression 
G(t;x; 7o ;#) = 0. 

Unconstrained in common we suppose that 

G(t; x; lo] 6) = G(t; x; 0) = 0. 



Such away, the system (13) which has an initial conditions y(t; x; 7) 



7=0 



and it has an uniqueness solving. This system we transform to the form: 



gy(s_7) 

d'y 



det 



ei e 2 

frx(t;y(-;7)) <9w(t; y(-; 7)) 

V3i(*;y(-;7)) ¥>32(*;y(-;7)) 

^i(*;y(-;7)) ¥ J n2(*;y(-;7)) 



gM(t;y(-;7)) 
V3n(*;y(-;7)) 
Vnn(t; y(s 7)) 



(7) 



Consequently, the function G(i;x;7) is automorphic transform for the func- 
tion u(t;x) and it has the next representation: G(t;x;7) = y(t;x;7) — x, 
where y(£;x;7) is a solution of the differential system (7) with an initial 

condition y(t; x; 7) 



x. 



7=0 



2 Construction of the program control 

On the analogy of the article [7] let us introduce the following definition. 



Definition 2 By the program motion of the stochastic system 

dx(t)= P(t;x(t)) +Q(t;x(t)) -s(t;x(t)) dt+ 
+B(t;x(t))dw(t) + J u(dt;dj)G{t;x(t);j), 



(18) 
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where w(£) is the m- dimensional Wiener process; u(t; A 7 ) - is homogeneous 
on t non centered Poisson measure, is meant the solution x(t; x G , s; u) which 
enables one under some (program) control s(i;x) and for all t to remain 

with the probability 1 on the given integral manifold u(t\ x(i; x )J = it(0; x Q ), 

which is the first integral of the Eq. (18) under the given initial conditions 



x(t; x ) 



t=o 



By this means we can construct the program control with the proba- 
bility equaled to 1 for dynamic systems which are subjected to the Wiener 
perturbations and the Poisson jumps. 

Theorem 2 The control s(t;x) allowing system (18) which is subjected to 
the Wiener perturbations and the Poisson jumps always remain on the dy- 
namically structured integral manifold u(t] x(t; x D ); coj = w(0;x o ) is deter- 
mined as the solution of the system of linear equation, which consists of the 
Eq. (18) and the Eq. (3). The coefficients of the second equation and respec- 
tively coefficient of the first equation of this system are determined by the 
Theorem 1. In addition we completely define the response to random action. 

Example 1 Needed is to determine the control s(t; x) and the response to 
random action enabling to dynamical system 



(x^t) + x 2 (t) + e~* + s 1 (t;x(t)))dt+ 



dxi(t) 
dt 

+b 1 (t;x(t))dw(t) + / 0i(t;x(t);7)i/(dfc;d7), 

JR( 7 ) 



(19) 

^i(t)x 2 (t) + e~ 2t + s 2 (t;x(t)))dt+ 
+b 2 (t; x(t))dw(t) + / g 2 (t; x(t); 7 )i/(dt; d 7 ), 

to remain on the given integral manifold w(t;x(t)) = x 2 (t)e~ 2xi ^ with prob- 
ability one. 

Solution. According to the Theorem 2 we construct the GSDE system, 
which has a function u(t; x) = x 2 e~ 2xi as the first integral. 

We first find an automorphic transformation of function u(t; x) = x 2 e~ 2xi . 
The corresponding partial derivatives of function u(t; y) = y 2 e~ 2yi are given 
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by 



Then we have 



dyi 



oy 2 



<9y(t;x; 7 ) 

<9 7 



Solution of this system is: 



d?/i(t;x;7) \ 
<9 7 

dy 2 (t;x; 7 ) 
d 7 / 



-22/1 



(20) 



y 1 (i;x; 7 ) = -ln(27 + 2C 1 (x)) 
</ 2 (*;x; 7 ) = C 2 (x) ( 7 + d(x)). 

With regard to initial conditions (according to Theorem 1, point 3), i. e. 

y(*;x;7) 



x, we obtain: 



7=0 



Ci(x) = -e 



2xi 



C 2 (x) = 2x 2 e 



-2xi 



And thus we have: 



1 



|/ 1 (t;x; T ) = -ln (2 7 + e 2 ^), 
y 2 (t; x; 7) = 2x 2 'je~ 2xi + x 2 . 

Consequently, automorphic transformation g(-) = (gi(-), of the func- 

tion x) = x 2 e~ 2:Cl is: 



#i(t;x; 7 ) = iln (2 7 + e 2 * 1 ) - x x , 
2 (*;x; 7 ) = 2x 27 e" 2:ri . 



(21) 



According to the statement 1 of the Theorem 1 we construct a matrix B (in 
this case it is a column as far as w(t) is one-dimensional Wiener process): 

B(-) = B(t;^) = q 00 {e' 2 *\ 2x 2 e~ 2 ^y , 

where q OQ = q 00 (t; x). Let us define the second expression in Eq. (5): 



"&B(t;x)~ 




<9x 


= Qoo ( 



-2e 



-2xi 



4x 2 e" 2xi 2e" 2:El ' ' 
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dB{t;x) 
<9x 



B(t;x) = q. 



-Ae~ 2xi 




-4q 2 00 e-^ 




det H{t; x) = det I -2x 2 e" 2xi e~ 2xi 



By Eq. (6) let's construct the matrix H(t; x) and their determinant 

e ei e 2 


/l /2 /3 

= e (-2/ 3 x 2 e- 2 ^ - / 2 e- 2 ^) + ei(/ ie -^) + 

+e 2 (2/ ia ; 2 e- 2 ^), 

where = /;(£; x), % = 1, 2, 3. 

In result the components of vector A = A(t; x) are given by 

h 



a i 



■ 4- 2a 2 e~ Axi 

/ 2 + 2/ S *2 + '"" ^ 



«2 



2/i* 2 



/ 2 + 2/ 3 * 2 



Hence, the desired GSDE system is: 

fx 



dxi(t) 



+q 00 e 



f2 + 2f 3 X 2 (t) 

(\ 



+ W 00 e 



2 -4an(i) 
oo 



- 2x ^dw{t) + / f- In (2 7 + e 2xiW ) - x 1 (t))i^(dt; dry) 



<ix 2 (t) 



2/1*2 (*) 



/2 + 2/ 3 * 2 (t)J 

+g o2* 2 (t)e- 2xiW rfw;(t) + / {2x 2 {t) 1 e~ 2x ^)v{dt; dry) 



JR(~f) 

Then we must solve the system of linear equation 

/l i o 7 2 --4*1 (t) 



+ * 2 (£) + e _t + s x (*; x(t)) 



*i(t)* 2 (0 + e~ 2t + s 2 (t;x(t)) 



/ 2 + 2/ 3 *2 (*) 

2/i*2(0 



+ 2<4e" 



Hence, program control is of the form of 

fi 



sx(t;x(t)) 



/ 2 + 2/3*2(0' 
+ 2g 2 e- 4a;i -x 1 (t)-x 2 (t)-e- 



/ 2 + 2/ 3 * 2 (t) 
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where f t = /j(t;x(t)), i = 1,2,3, q OQ = q 00 (t;x(t)) and f 2 + 2f 3 x 2 (t) ^ 0. The 
response to the Wiener action is defined as 

6i(t; x(t)) = gooe- 21 ' 1 ^), 6 2 (t; x(t)) = g 00 2x 2 (t)e- 2:El « 

Components of the response to the Poisson jumps are defined as: 

^(t;x(t); 7 ) = iln (2 7 + c **<*>) _ ^(t), 
p 2 (t;x(t); 7 ) =2x 2 (t) 7 e- 2 ^. 

Therefore, the set of vectors of the program controls and their corre- 
sponding responses to a random action enabling the dynamical system under 
consideration to remain on the given time-variable manifold are determined. 

The choice of functions /j(£;x(£)), % = 1,2,3 and g 00 (i;x(t)) allows to 
construct the program control on account of some conditions. For example 
we take into account a simulation utility. 

Remark 2 It is worthy to note that the manifold can be defined by a set of 
functions [6]. 

Conclusion 

Application of the first integral stochastic theory allows to construct the 
program control for dynamic system with probability equaled to 1 of there 
are strong random perturbations, caused by Wiener and Poisson processes. 
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